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We present results from a covariance study for the proposed Laser Astrometric Test of Relativity 
(LATOR) mission. This mission would send two laser-transmitter spacecraft behind the Sun and 
measure the relative gravitational light bending of their signals using a hundred-meter-baseline 
optical interferometer to be constructed on the International Space Station. We assume that each 
spacecraft is equipped with a < 1.9 x 10 _13 ms 2 Hz _1//2 drag-free system and assume approximately 
one year of data. We conclude that the observations allow a simultaneous determination of the orbit 
parameters of the spacecraft and of the Parametrized Post-Newtonian (PPN) parameter 7 with an 
uncertainty of 2.7 x 10 -9 . We also find a 5 x 10~ 9 determination of the solar quadrupole moment, 
J2, as well as the first measurement of the second-order post-PPN parameter S to an accuracy of 
about 10" 3 . 



I. INTRODUCTION 



In this paper, we present results of a covariance study for the proposed Laser Astrometric Test of Relativity 
(LATOR) mission (see Turyshev et al. [|). The goal of the LATOR Mission is to effect a precise measurement of 
the gravitational bending of light for laser signals passing near the Sun and thereby determine several Parametrized 
Post-Newtonian (PPN) parameters of relativistic gravity with an accuracy many orders of magnitude better than 
\^ ' present best estimates. 

The LATOR instrument consists of two Sun-orbiting laser transmitter spacecraft in eccentric ecliptic orbits with a 
3/2 period resonance with the Earth's orbit. The two spacecraft are separated by about 1 degree, as seen from the 
Earth. The angular position of each spacecraft is measured using a lOOm-baseline optical interferometer, proposed for 
■ the Earth-orbiting International Space Station (ISS). Figure^shows the LATOR mission geometry. The coherent laser 
signals from each spacecraft will allow the position of each to be determined with an accuracy of 0.1 picoradians (20 
CJ nanoarcseconds) , relative to the reference frame defined by the interferometer baseline. The simultaneous measurement 
^ ' of two spacecraft enables a relative angle to be measured, eliminating the need to tie the interferometer frame to inertial 
?H . space. As the lines-of-sight from the spacecraft to the Earth pass close to the Sun, the gravitational bending of light 
^P" will affect the inner spacecraft's apparent position more than it does the outer spacecraft, causing the inter-spacecraft 
angle, as seen from Earth, to differ from its value in the absence of gravitational effects. Measurements of this inter- 
spacecraft angle may thus detect the gravitational bending of light and verify gravitational theories that predict it. 
Many measurements of this angle will be made as, over the course of a year, the spacecraft will make three such passes 
behind the Sun, giving detailed information about how the bending of light changes depending on which part of the 
solar gravitational field the light passes through (see figure |5J). For more information on the design of LATOR, see 
Turyshev et al. @,GJ. 

One of the problems, of course, in measuring the gravitational deflection of the laser signals from LATOR is to 
know what the apparent angle between the two spacecraft would be without the gravitational bending present. Two 
methods have been proposed to provide this knowledge. One is to have a microwave (or optical) signal sent between the 
two spacecraft, determining their separation to an accuracy of ~ 1 cm. The other is to equip each LATOR spacecraft 
with a drag-free system, so that the orbit of each is very close to a purely gravitational and noise-free trajectory. The 
spacecraft-to-spacecraft distance can then be inferred by accurate tracking and orbit-modeling over the time of the 
mission. As we will discuss in section IIVI the spacecraft-to-spacecraft link leaves one important problem unsolved, 
that of determining the impact parameter. However, as we show by means of this covariance study, the drag-free 
system allows orbit modeling using the primary LATOR data to simultaneously solve the spacecraft-to-spacecraft 
distance and the impact parameter problems. 
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FIG. 1: Diagram of LATOR geometry: The transmitting spacecraft are on far side of the Sun from Earth; the receiving 
spacecraft are in low Earth orbit aboard the ISS 



II. LATOR SCIENCE GOALS 



As was shown in the early 1970's, most theories of gravity are metric theories in which the dynamics of neutral matter 
are determined by a metric, which is itself determined by the neighboring mass distribution. Will and Nordtvedt 
developed an expansion of the metric that has a tunable set of numerical parameters which allow it to describe 
a very broad range of metric theories of gravity. This Parameterized Post-Newtonian (PPN) expansion is useful for 
experimental tests of gravity because the task of fitting the model to the data is carried out by simply adjusting the 
numerical values of these tunable parameters so that the model best fits the data. For the LATOR experiment, the 
version of the PPN metric required is: 

GM . R 2 3cos 2 (9-L nn G 2 M 2 
500 = -1 + 2— (1-J 2 -^ — (!) 

r „ o GM u r i? 2 3cos 2 0-l. 3 G 2 M 2 
9lm = Si m {l + 2 1 —[1-J 2 - T ] + _£__), (2) 

where M is the mass of the Sun, R is the solar radius, r and 9 are polar coordinates, J2 is the solar quadrupole 
moment, and 7, f3 and S are PPN parameters. While the full PPN metric contains other terms as well, this version 
contains the terms that contribute significantly to the bending of light that is observable by LATOR. There is a 
frame-dragging term that scales with the solar rotational angular momentum, J z , and which should have an effect 
similar in size to that of other terms observable by LATOR (see Richter and Matzner [f| ) . This term is not included 
in our model, but should not be difficult to add at a later date. Table [I] contains a description of the effects of 
each parameter in the above metric, along with its current best estimate and theoretical value in Einstein's General 
Relativity (GR). 

The accuracy of an interferometer angular measurement of gravitational deflection is related to the accuracy of the 
measurement of the differential time-of-arrival of a phase front in the two detectors of the interferometer. The relation 
is 

L 

where A9 is the accuracy of the angular measurement in radians, At is the accuracy of the timing measurement 
in seconds, and L is the interferometer baseline in meters. LATOR's 20 nanoarcsecond (or 10 -13 radian) angular 
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FIG. 2: Plot of impact parameter vs. time for orbits of spacecraft 1 (crosses) and spacecraft 2 (diamonds), with 1 day sampling 
interval. Time axis is in days, vertical (impact parameter) axis is in units of Rsun- 



Parameter 


Significance 


Current best estimate 


7 


Space-time curvature per unit 
rest mass. 


1 + (2.1 ±2.3) x 10 _5 |3. Equal to 
one in GR 


P 


Space curvature per unit gravita- 
tional self-energy 


1 + (0.9±1.1) x 10" 4 Q. Equal to 
one in GR 


8 


Time curvature per unit gravita- 
tional self-energy 


Unmeasured. Equal to one in GR 


■h 


Solar quadrupole moment param- 
eter. 


Never measured directly. Esti- 
mated from solar models to be 
about 10" 7 . 



TABLE I: Significance and best estimates for parameters in PPN metric 



accuracy thus corresponds to 3.33 x 10~ 20 s (or about 20 /^.cycles) accuracy in timing measurement at each element of 
the interferometer. To second order in GM/r, the time delay derived from the metric of eqs. Qandl^lmay be written 
(see Richter and Matzner 

ttr x R -x T GM R 1, 2 {GM) 2 R 2 x T 2 x R x T ~y\ 

—5- = — 5 — I 1 b h n 1 + 7 — 52 rU 1 ) + 2 ) 

R R R tt 2 R tt TbTt 



Sun 



FIG. 3: Diagram of coordinate system used for equation [3] 




(arctan — — — arctan — — ) + 

Vt Ut Vt 



n r R x T x R 1 

(1 )( — + 1) + -(1+7)J 2 

r T r T r R 2 



GM R 2 



(3) 




Here, it is assumed that the x axis is tangent to the photon trajectory at the transmitter, the origin is at the position 
of the source mass (the Sun), and that the Sun, transmitter and receiver are in the z = plane. R is the radius of 
the Sun, x R and r R are the receiver x-coordinate and distance from the origin, xt and yx are the transmitter x- and 
y-coordinates, and rr is the transmitter distance from the origin (see figure [3J). 

The sizes of each term in eq. |3] for photons travelling 2 AU and just grazing the limb of the Sun are given in 
table In] also taken from Richter and Matzner Q . The sizes of these terms should be compared with the timing 
accuracy expected for the LATOR interferometer measurement, remembering that the actual differential times-of- 
flight measurements for the two LATOR signals are a fraction of the displayed terms in table ITU 

The leading term, after the 2000-second geometrical time-of-flight term, is the term proportional to 1 + 7, of 
size 2.4 x 10 -4 , with its second-order correction of size 3.6 x 10 -8 . This is the effect whose measurement will give 
PPN parameter 7 from the LATOR measurements. The next largest term, with size 4.9 x 10 -10 enables the solar 
quadrupole moment, J2, to be measured, while the term of size 2.5 x 10~ 10 will give a linear combination of 7, j3, and 
the post-post-Newtonian term 5. 

The fundamental current incompatibility of General Relativity and Quantum Mechanics suggests that the PPN 
parameters should deviate from their values in GR at some level. The size of the deviation of the PPN parameters 
from their value in GR is theory dependent. Damour and Nordtvedt 0,0 suggest that for some scalar-tensor theories 
of gravity the deviation in the PPN parameter 7 occurs at the level 1 — 7 ~ 10~ 7 , significantly smaller than the 
uncertainty in the current best estimate for 7, at ~ 10~ 5 (see Bertotti et al. 6]). The primary goal of the LATOR 
mission is measurement of the PPN parameter 7 to an accuracy of 10 -9 . Secondarily, LATOR will provide the first 
measurement of a post-post-Newtonian parameter, S, and a measurement of the solar quadrupole moment parameter 
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Term: 


Size of term in seconds: 


4£ 


2000 


2(l+ 7 )^ln(4(£) 2 ) 


2.4 x 10~ 4 




2.0 x 10" 5 


4(1+7)W) 2 £ 


3.6 x 10~ 8 


2(1 + 7)J 2 — 


4.9 x 10" 10 


4[2(l+ 7 )-/3+!<5](^) 2 arctan(^) 


2.5 x 10~ 10 


2(l+7)(^) 2 ^ln(4(£) 2 ) 


2.4 x 10~ 12 


2(2/3 -l) ( eM) 2 -£ 


9.7 x 10" 14 



TABLE II: Size of terms in time-of-fiight expression for 2 AU light path which grazes the limb of the Sun. PPN parameters 
are set to their values in GR, and J2 and J z are set to their values in Dicke's model of the Sun. 



Parameter 


Description 


ai, 02 


Semi-major axes for orbits of transmitting spacecraft 1 and 2. 


ei, e 2 


Eccentricities for orbits of transmitting spacecraft 1 and 2. 


001, <t>02 


Initial phase for orbits of transmitting spacecraft. 


U)l, W2 


Angle of semi-major axes of transmitter orbits with respect to x axis. 


0x1, a>±2 


Constant acceleration in drag- free system, perpendicular to transmitter-receiver line-of-sight, 
for transmitting spacecraft. 


a ||l, a ||2 


Constant acceleration in drag-free system, parallel to transmitter-receiver line-of-sight, for 
transmitting spacecraft. 


7 


Space-time curvature per unit rest mass. 


8 


Time curvature per unit gravitational self energy. 


J2 


Solar quadrupole moment parameter. 



TABLE III: Parameters to be solved for in covariance study 



J2, which will provide information about the solar interior. 

III. COVARIANCE STUDY 

We model the LATOR instrument in the solar system and use a linear least squares covariance analysis to solve 
simultaneously for the uncertainties in the PPN parameters and orbit parameters. A Newtonian model is used for 
the transmitter orbits, and time delays for the interferometer legs are modeled using eq. |21 A priori uncertainties for 
the orbit parameters are assumed, consistent with typical uncertainties in spacecraft radio tracking. The transmitter 
coordinates are evaluated as functions of the time and of the orbit parameters for each spacecraft, while receiver 
coordinates are evaluated as functions of time using a simple circular heliocentric orbit. A justification for this simple 
model of the transmitter and receiver orbits is given in appendix El The interferometer is assumed to be inertially 
oriented. We model two constant body-fixed accelerations for each transmitter, one along the Earth line-of-sight and 
one perpendicular to it. For each assumed data point, the time-of- flight from each spacecraft to each antenna is 
calculated as a function of the time of observation and of the parameters of the model. The parameters include the 
constants of motion for the transmitter orbit (one set for each transmitter), the DC accelerations on each spacecraft 
(one set per transmitter) and the PPN parameters. These parameters are listed in table lllll The goal of the covariance 
analysis is to estimate the uncertainties with which these parameters will be determined in the final LATOR mission 
data analysis. 

Notably absent from the list of parameters is the parameter 0. This is because the only significant term in the 
time-of- flight equation which contains (see table [HJ is 

4[2(l + 7 )-/?+^](^) 2 arctan(^) 

in which 0, S and 7 are all completely correlated. Other terms give 7, but the only other term to contain is smaller, 
meaning that would be less well determined than the combination + 3/45 and that the incertainty in both would 
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be that of j3 alone. However, since j3 is already known from other solar system experiments (Williams et al. Q) 
with an accuracy (see table better than the best LATOR can manage for this term, our analysis assumes that (3 is 
already known. The term 

_ 3 , , (jt J\/j[ o , r , 

4[2(1 ^](—f arctan(-) 

is therefore used to solve for S alone. 

At each time sampled, the LATOR signal is formed from 

S = [Ti2 - Tn] - [t 2 2 - T2l] (4) 

where each t is the time-of-flight for a particular transmitter-receiver pair; so, for example, t\% would be the time- 
of-flight between transmitter 1 and receiver 2. At each time, the contribution from the current time step to the 
information matrix is formed by calculating the derivatives of S with respect to the parameters of table ITTll Referring 
to the parameter set collectively as x, the contribution to the information matrix at each time U is given by 

1 dS(U,x) dS{U,x) ... 
a ' k{U) = W~tei dxV~ (5) 

where a is the uncertainty in LATOR's measurement of each time delay. In our analysis, a was set to 3 x 10~ 20 s, 
corresponding to the conservative 10 picometer interferometer uncertainty estimate found in LATOR literature (see, 
for instance, Turyshev et al. 0). 

The derivatives in eq. were calculated using analytic derivatives for PPN parameters, finite differencing for 
orbit constants of motion, and finite differencing with a perturbative Enke integration for the DC accelerations. 
Calculations leading up to the time-of-flight calculation were performed in double precision. Round-off error was a 
serious problem for derivatives with respect to orbit parameters, so we used Ridders' extrapolative finite differencing 
algorithm (from Press et al. 01 chapter 5) and performed all calculations after the time-of-flight calculation in the 
GNU Multiple Precision Arithmetic Library (CMP) in order to keep sufficient precision to solve for uncertainties in 
all of the variables. 

The terms in the information matrix, ajk(U), were then summed over a period of about 1 year, during which the 
transmitting spacecraft pass behind the Sun three times as viewed from Earth (see figure [2J . We used a 1-minute 
sampling cadence, but the number of data points was further reduced because the transmitting spacecraft are only 
visible from Earth for half of the ISS hour-and-a-half orbit, and it may take tens of minutes per orbit for the receiving 
interferometer to acquire lock. This was simulated by adding to the information matrix for only 10 minutes out of a 
90 minute period. The resulting integration has about 50000 time steps. The total information matrix is 

^ 1 dS(U,x)dS(U,x) 

^ = 1.^—^ — (8) 

The information matrix was inverted using a Singular Value Decomposition (SVD), and the diagonal elements of 
the inverse gave the standard variations for the corresponding parameters: 

(no sum) (7) 

IV. THE IMPACT PARAMETER PROBLEM 

Since the leading 7 term in the time-of-flight formula is logarithmic, its uncertainty scales as 1/r. We will so far 
anticipate the results of our analysis to say that we expect 7 to be determined to an accuracy 3 x 10~ 9 . If this is 
to be accomplished, then the size of the coefficient of 7 must be known correctly to the same accuracy. For signals 
passing close to the Sun, the 1 + 7 term in eq. OH may be approximated as 

MM, GM R . , x R + r R .„ .GM, .„ . 

(l + 7)-p-(l"-p In , ~ 1 + 7 — M4r fl r T /d 2 (8) 

K K tt xt + tt H. 

where d is the impact parameter. Thus, if the coefficient of 7 is to be known to 3 x 10 9 , the impact parameter must 
be known to a part in 10 9 . For paths passing near the solar limb, the impact parameter is ~ 10 9 m, meaning that it 
must be known to within a meter. How is this to be done? 
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First, we point out that the spacecraft-to-spacecraft link that could be used to give the base of the narrow triangle 
linking the two spacecraft and the interferometer provides no information on how that triangle is oriented relative 
to the Sun. Observations of the solar limb using the interferometer cannot possibly be made accurate to 1 meter, 
and, even then, the distance from the solar limb to the center of the Sun (for purposes of calculating the solar 
potential) cannot possibly be modeled. One could conceive of doing orbit determination using some typical Doppler 
or ranging tracking of the spacecraft or even of using the precise relative angular measurements from the interferometer 
themselves as a new navigation data type. However, experience with orbit determination for free-flying interplanetary 
spacecraft shows that non-gravitational forces do not allow anything like 1-meter accuracy over periods of months. 
If the impact parameter problem is to be solved, then non-gravitational forces must be eliminated using a drag-free 
system on board each spacecraft. 

In our covariance study, we have assumed that the LATOR spacecraft follow purely gravitational orbits (except for 
a constant body- fixed acceleration that is solved for in the covariance study). We will find in our study that not only 
is the impact parameter problem solved by the subsequent orbit determination, but that the spacecraft-to-spacecraft 
distance problem is likewise resolved. This is to say that, as we simultaneously solve for orbit parameters for both 
spacecraft along with the PPN parameters, any uncertainties in the inertial positions of the spacecraft are correctly 
taken into account through the correlations in the covariance matrix. Thus, whatever final uncertainty we derive 
for the PPN parameters will be the correct uncertainty, given the uncertainties in the data-derived orbits of the 
spacecraft. 

To determine how good a drag-free system needs to be to enable a pure gravitational trajectory to be a good 
approximation, let us consider the following: 

1. A 10 _13 -radian angular accuracy at a distance of 3 x 10 11 m corresponds to a position accuracy of 3 cm. 

2. The position noise for a spacecraft suffering acceleration noise with noise spectral density S a is x lms — 
\J 5' a /(167r 4 / 4 ), at a Fourier frequency /. 

3. Therefore, if we require 3 cm accuracy at a frequency 1/month (4 x 10~ 7 Hz), we derive a drag-free requirement 
of S a = 3.6 x 1(T 26 m 2 s -4 Hz -1 , or a root spectral density y/S2 = 1.9 x 1CT 13 ms^Hz" 1 / 2 . 

The currently planned LISA reference sensors are expected to have an accuracy of 3 x 10 ms" - 2 Hz _1/2 . Below 
10 Hz, a 1// increase is expected due to thermal drift resulting from the passive thermal isolation of the reference 
sensors. At 3 x 10 _8 Hz, this would produce acceleration noise of 10~ n ms~ 2 Hz -1 / 2 , so a mK active thermal control 
(not a difficult technology) would be needed to keep thermal noise sources below the 2 x 10~ 13 ms -2 Hz -1 / 2 we have 
assumed for the LATOR accelerometer. Voltage drifts in the measurement electronics do not affect the capacitance 
bridge used to read the accelerometer, since it is nearly a null measurement. Proof-mass charge control is absolute and 
independent of long-period effects. Thermal drift of the non-controlled parts of the spacecraft, along with mechanical 
relaxation of the spacecraft structure, will produce long-term self-gravity changes, but we have calculated these to be 
well out of the way of the level of performance we have assumed, as long as there is a level of spacecraft mechanical 
and thermal engineering comparable to that already required for LISA. Most other noise sources are white and LISA 
performance would keep them below the 10~ 15 ms -2 Hz -1 / 2 LISA requirement. 

V. RESULTS AND CONCLUSIONS 

The predicted accuracies of orbit parameters and spacecraft DC accelerations from one year of LATOR data are 
shown in tabic |IV| The 60 -cm uncertainty in the two semi-major axes is primarily the result of strong correlation 
between these parameters. (The uncertainty for a% or ai would have been about 60 microns if one of these two 
parameters were being solved for alone.) 

The results for the PPN parameters are shown in table ]V\ The uncertainty in 7 will represent a four order-of- 
magnitude improvement over the current Cassini results and will, for the first time, approach the level where heuristic 
string-theory arguments predict that a deviation from GR might be expected. The Parametrized Post-Post-Newtonian 
(PPPN) parameter S will be measured for the first time. This is important because there are several theories whose 
PPN parameters are identical with GR, but which differ from GR in their predictions for the PPPN S term. Finally, 
we predict a 6% determination of the solar quadrupolc moment, relative to its expected value of 10 -7 . This is a 
direct, model-independent value and may be used as a touchstone value for various models of the solar interior. 

We conclude that LATOR will represent a significant step forward in the experimental tests of relativistic theories 
of gravity and will be able to satisfy the goals that were suggested in the initial project proposals. We find that, with 
a 1.9 x 10~ 13 ms -2 Hz -1 / 2 drag-free system covering a band down to nearly DC, it will be possible to simultaneously 
extract the PPN parameters and all relevant orbit parameters from the interferometer data alone. This seems to be 
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Parameter 


LATOR predicted uncertainty 


ai, d2 


1.4 m 


ei, e2 


2 x 10" 11 


001, <f>02 


1.3 x 1CT 10 


Wl, CJ2 


4 x KT 11 


a±i,a±2 


4 x 10~ 14 m/s 2 


a ||l, a ||2 


3 x 1(T 14 m/s 2 



TABLE IV: Covariance study estimates for LATOR determination of uncertainties in transmitter constants of the motion and 
DC accelerations 



Parameter 


Significance 


Current best estimate 


LATOR predicted uncer- 
tainty 


7 


Space-time curvature per unit rest 
mass. 


1 + (2.1 ± 2.3) x lCr 5 Q. Equal to 
one in GR 


2.7 x 10 -9 


P 


Space curvature per unit gravita- 
tional self-energy 


1 + (0.9 ± 1.1) x 1O" 4 0- Equal to 
one in GR 


Not solved for. 


8 


Time curvature per unit gravita- 
tional self-energy 


Unmeasured. Equal to one in GR 


1 x 10" 3 


J 2 


Solar quadrupole moment parame- 
ter. 


Umeasured. Estimate of about 
10" 7 inferred from solar models. 


5 x 10~ 9 



TABLE V: Predicted uncertainties in PPN parameters with LATOR and comparison with current best estimates. 



the only good solution to the impact parameter problem, and we strongly suggest that the engineering studies for 
LATOR consider this hardware option. 

APPENDIX A: JUSTIFICATION OF ORBIT MODELLING 

The assumed accuracy of each interferometer measurement is 10 -13 radians, which, at a distance of 2AU (3 x 10 11 m) 
corresponds to 3 cm along-track position accuracy. As a comparison, let us consider the recent range data to various 
Mars spacecraft, the forces on Mars being purely gravitational. The data for the last decade have post-fit residuals of 
1-2 m. The parameters that enter the solution include spacecraft locations and Mars physical ephemeris for landers, 
orbital elements and Mars gravity field for orbiters, corrections to Earth tracking station locations, solar corona and 
troposphere delays, orbits and masses of the planets, and the masses of several asteroids. This combination of ~ 200 
parameters in comparison with our 4 orbital elements would seem to suggest that we have oversimplified the problem. 
However, the LATOR data and model are immune to the effects of most of these elements of the model. 

The LATOR interferometers are in orbit above the troposphere. The laser frequencies are high enough to avoid 
significant plasma propagation delays. There is no local disturbing planet to complicate the spacecraft motion. The 
interferometer makes an angular position measurement, so errors in the interferometer position cancel to first order. 
Errors in the interferometer orientation likewise cancel when the differential position of the two spacecraft is calculated. 
Finally, the long-period solar system model uncertainties in the position of each spacecraft are significantly reduced 
when the differenced spacecraft positions are calculated. For an along-track sinusoidal error for each spacecraft of size 
s and period T (in days) , the size of the error in the differenced position is given by 

As = '1^1/Ilsx (Al) 
vT 

where v is the velocity of the spacecraft in degrees per day and SX is the difference in heliocentric longitude between 
the two spacecraft, in degrees. Since v *~ 1 and 5X ~ 2°, and since the solar system model errors are of order 
~years, the error in the differential position of the two spacecraft will be reduced by a factor ^ (which is ~ 0.01 
or less) relative to the single-spacecraft position error. Thus a spacecraft error of 1-2 m becomes an error of 1-2 cm 
in the LATOR observable, slightly below the 3 cm error we assume in our covariance study. It is therefore not an 
unreasonable assumption to treat the solar system parameters as known to this level, by virtue of the present and 
ongoing accurate Mars ranging data. 
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